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! Q \ The general solutions of many three-dimensional Lotka-Volterra systems, previously 
^— i ■ known to be at least partially integrable, are constructed with the aid of special functions. 
I/-) \ Examples include certain ABC and May-Leonard systems. The special functions used 
are incomplete beta and elliptic functions. In some cases the solution is parametric, 
with the independent and dependent variables expressed as functions of a 'new time' 
variable. This auxiliary variable satisfies a nonlinear third-order differential equation 
of a generalized Schwarzian type, and results of Carton-LeBrun on such equations are 
exploited. Several difficult Lotka-Volterra systems are successfully integrated in terms of 
Painleve transcendents. An appendix on incomplete beta functions is included. 
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1. Introduction 

(a) Background 

Models of Lotka-Volterra type occur frequently in the physical and engineering 
sciences, as well as the biological. In any such model there is a (i-dimensional state 
vector x = (x\, . . . , Xd), a function of time r, which satisfies an autonomous system 
of ordinary differential equations (ODE's) of the form 

Xi = Xi ^a i0 + ^ aijXj^j , i = l,...,d. (1.1) 

The overdot signifies d/dr. By letting xq = 1 and ao« = one can optionally rewrite 
this as Xi = Xi X/?=o a ij x ji * = 0, . . . ,d. There is a paucity of closed- form solutions: 
when d ^ 3 and even when d = 2, the system (jl.ip is usually integrated numerically 
rather than symbolically. 

In most applications Xi ^ for all i, though it may be useful to allow the Xi 
and eve n r to be complex. Such models were first introduced in population 
ecology ( Kotll200ll : lMavll200ll ). where x lies in the non-negative orthant, and the 



terminology reflects this: a^o is the growth rate of the i'th 'species,' and depending 
on the signs of the elements of the interaction matrix A = (aij)fj =1 , one speaks 
of predation, competition or mutualism. If a^o > and an < 0, in the absence of 
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the other d — 1 species the i'th species will grow logistically, but the behaviour of 
x = x(t) in the interior of the orthant may be much more complicated. 

Modelling by small-(i Lotka-Volterra systems has occurred many times in 
physics. I n laser phy sics this ranges from the initial treatment of multimode 
coup ling ( Lamb|[l964 ) to Raman amplification in an optical fibre (ICastella et al. 



I2008h . where r is the distance along the fibre. Other physical applications include 
the modelling o f the (integrable) interaction of Langmuir waves in plasmas 
( Qin et al. 201 ll ) and more chaotic phenomena: the Ku eppers-Lortz instability in 
a rotating fluid f|Busselll98.sl : ISan Miguel fc Tor all 120001 ) . for which d = 3, and the 
Boltzman n dynamics of d rarefied, spatially homogeneous gases in a background 
medium (Jenks 1969; iLupini fc Spiga 1988 ). Several integrable Lotka-Volterra 
systems with fairly small d, of mathematical interest, have been obtained as 
spatial d iscretizations of the Korteweg-de Vries eq uation, or as generalizations 
of same (jBogovavlenskv 1988 : Bogovavlenskii 2008 ). 

In population biology, small-c? Lotka-Volterra systems include the classical 
d = 2 predator-prey model of Volterra and Lotka, well-known though structurally 
unstable; the famous May-Leonard model of d = 3 cycli cally competing specie s 
with equal growth rates, which has a stable limit cycle (IMav &: Leonard! Il975[ ): 
and genera lizations that model the asymmetric com petition of d species in a 
chemostat ( Wolkowicz 2006 : Ajbar &: A lhumaizil l2012l ). There is also a literature 
on evolutionarily stable strategies in game dynamics, modelled by diff e rentia l 
systems equivalent to Lotka-Volterra ones (jHofbauer &: Sigmundl Il988l . |2003| ) . 
Lotka-Volterra systems with d^3 may display chaotic behaviour, a fact first 
noted in ecological modelling ( Klebanoff &: Hastings! T994I : Vano et al. 2006 ). 

Many nonquadratic nonlinear differential systems ca n be transformed 
to (11.11) by chang es of variable (IPeschel Sz Mendel Il986l : iBrenig Gorielv 

Complex chemical 



iges ( 

19891 : iGouzi Il992l : iHernandez-Bermeio Fairenl 



reactions with mass-action kinetics provide examples. For special stoichiometries, 
spatially homogene o us species concentrati ons evolve in a Lotka-Volterra way 
fjErdi fc Tothl Il989l : iMurza fc Tercel l2010l ). Furthermore, mass-action systems 
with nonquadratic polynomial nonlinearities, which can be obtained from (formal) 
chemical reactions in more than one way, can often be transformed to the Lotka- 
Volterra form. 

Several elementary changes of variable deserve mention, for any cl-dimensional 
system of the form (jl.ip . First, by scaling r one can scale the column vector of 
growth rates and the interaction matrix A. Additionally, by scaling x\,...,Xd 
independently one can scale independently the d columns of A. Thus the 
parameter space has dimensionality not d 2 + d but d 2 — 1. 

If the growth rates a\o, ... ,0,40 are equal (to a*o, say), then in terms of a 
transformed time r 1 :=e a *° T the transformed variables x\ := e~ a *° T Xj will satisf y 
a system of the form (jl.ip but with zero growth rates. Furthermore ( Planklfl999l ). 
if x\, . . . , Xd satisfy (jl.ip with zero growth rates then the ratios (x\, . . . , x^-i) '■= 
(xi/xj, . . . ,Xrf_i/x^) will satisfy a (d — l)-dimensional system resembling (jl.ip . 
but with growth rates ajo := o-id ~ a dd aiK1 interactions dij := aij — a^j ■ Thus 
Lotka-Volterra systems with equal growth rates, for which the parameter space 
dimensionality is effectively d 2 — d, should be peculiarly amenable to analysis; 
though even they have usually been studied not symbolically but numerically. 
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The object of this paper is the construction of general solutions x = x(r) of 
systems of the form (jl.ip . with the aid of special functions. The focus is on d = 3 
systems, though ancillary results on integrable d=2 systems are obtained first. 
The special functions employed are the incomplete beta function and its inverse, 
which sometimes reduces to an elementary or elliptic function. For many d = 3 
systems known to be at least partially integrable, general solutions are constructed 
for the first time, with unexpected ease. By default, each of x, r is expressed as a 
function of an auxiliary variable t, a 'new time.' A minor but interesting example 
is the fully (rather than merely cyclically) symmetric May-Leonard model. 

Most of the integrated systems have equal growth rates, but not all. Several 
systems with unequal ones which have the Painleve property, that each solution 
x = x(t) can be extended analytically to a one-valued function on the complex 
r-plane, are integrated with the aid of Painleve transcendents. 

Many of the interesting examples of general solut i ons x = x(r) that are 
presented are made possible by results of Carton-LeBrun ( 19691 ). She classified all 



nonlinear third-order ODE's of a certain type that have the Painleve property; and 
the ODE satisfied by the new time t = t(r) used here was fortuitously included. 
She also integrated many such ODE's with the aid of elliptic functions. The 
examples in §§[3] and [4] may serve to reawaken interest in her work. 

This paper includes an appendix on closed-form expressions for inverse 
incomplete beta functions, which should be of independent interest. 

(6) Previous results, reformulated 

The following is a summary of symbolic (non- numerical) results on small-d 
Lotka-Volterra systems, which provides the context for the new results. 

Lotka-Volterra solutions cannot always be expressed in closed form, i.e. 
in terms of elementary or 'known' functions, even if the behaviour of the solutions 
is nonchaotic; or so it is believed. An example is the original system of Volterra and 
Lotka, with d = 2 and no self-interactions (an = 022 =0). If the growth rates are 



equal it has a general solution in terms of elementary functions (|Varmalll977l ). If 
the rates are unequal it is nonetheless integrable (see below) and hence nonchaotic, 
but this does not imply the existence of a useful closed form for the general 
solution. The case of unequal growth rates when d = 2 has in fact been u sed as a 



test-bed for geometric numerical integration schemes (lHairer et al.ll2006f ). 

Owing to the difficulty of constructing explicit solutions of systems of the 
form (jl.ip , there has been much work on the simpler problem of finding constants 
of the motion, i.e., first integrals. Any such is a function I = I(x\, . . . ,xa), 
preferably one-valued and well-behaved, obeying 1=0 on any system trajectory. 
(By the overdot the total time derivative is meant.) If such an / exists, any 
trajectory is confined to a [d — l)-dimensional surface 1= const; which, e.g., 
facilitates the study of long-time behaviour. If there are d — m functionally 
independent first integrals then motion is confined to a m-dimensional surface, 
obtained as an intersection. When m=l the system is completely integrable (its 
general solution can be 'reduced to quadratures,' though not necessarily usefully). 
When m = 2 it is partially integrable and is at least nonchaotic. 

For d = 2, 3, restrictions on the Lotka-Volterra parameters {aij}f =1 j =0 that 
imply the existence of a single closed-form first integral have been investigated by 
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many authors. In most cases the resulting first integral is of the form 



;i.2) 



with each fk = fk{ x l: ■ ■ ■ i x d) a polynomial (usually though not always of 
degree 1 with constant ter m allowed). Result s along this line include those 
of Cairo and collab orators (ICairo fc Feixl Il992l : ICairo et all Il999l : ICairol 1200(11 : 
ICairo &: Llibre 20001). who in many cases used the classical technique of Darboux 
polynomials (lGorieryll200ll) : and results based on an integrating factor technique 
(jSaputra et al.M2010f ). When d = 3 and a single first integral of a form similar 



to (|1.2p exists, there has been some work on the construction of a second, 
functi onally independent first integral of a more complicated and less algebraic 
form ( Grammaticos et al. 19901 : Gorielv 1992 : Gao 2000l : Bustamante fe Hoiman 



2003). 



The Darboux polynomial (DP) concept is fundamental. A DP / = 
f{x\, . . . , Xd) is a polynomial for which / = Kf along any system trajectory, where 
the 'cofactor' K = K(x±, . . . , Xd) is necessarily also a polynomial, of degree 1 with 
constant term allowed. The algebraic surface / = is thus invariant under the 
Lotka-Volterra flow. If the cofactors of DP's fi, ■ ■ ■ f r are linearly dependent, 
some product of the form (j 1 . 2 [) will be a first integral: generically, a non-trivial one. 
Finding a first integral, including the needed parametric restrictions, is facilitated 
by the coordinate functions x\,...,Xd being DP's, associated to the invariant 
hyperplanes Xj = 0, i = 1, . . . ,d, and having respective cofactors Ki = X^=o a ij x j- 
To construct a first integral a single additional DP / will suffice, if one sufficiently 
restricts parameters to ensure dependence of the cofactors of x%, . . . , x^, /. 

Tables Q] and [2] are a distillation of the results of the preceding authors on 
degree-1 DP's / that yield first integrals of the form \xi\\ . . \xd\ ld In the 

two tables, i, j and i, j, k are arbitrary permutations of 1, 2 and 1,2,3, respectively; 
and / = is an invariant line, resp. plane. For both d = 2 and 3, the cases when 
^ d restrictions on the parameters yield a first integral of this form split int o three 
types: I, II and III. ( The numbering here d iffers from that of lCairo &: Feixl (jl992l ) 
but agrees with that of iHua et al.l (1996).) The primed cases can be viewed as 
projectively transformed versions of the unprimed ones. For instance, the d = 2 
case II- is related to II by (xo,Xi,Xj) = (xj,Xi,xo)/xj. 

The type-I cases are degenerate: the first integral / does not require any 



additional DP / and is merely of the form F|° 



In fact, in each type-I 



case the restrictions on the parameters, expressed in terms of the d + 1 column 
vectors aj := (a,ij)f =l , j = 0, . . . , d, imply that the matrix A = (a,ij)f =1 - =1 must be 
singular, i.e., det A = 0. Lotka-Volterra models with singular A (e.g., ones in which 
A is antisymmetric and d is odd) may have physical applicatio ns, though their 
relevance to pop ulation ecology has been strongly questioned (iMavnard Smith 



19741 : lMavll2001[ ). But they will not be considered further here. 

For d = 2, in case II (which is defined by a±Q = 020 = 0, i.e., by the growth rates 
of the two species being zero), the resulting first integral is easily seen to be 

r_|„ |a22(a2i-a.n) |„ (0.12-0.22) \( n „ \„ ( „ ,, \„ |A /-, q\ 

1—\X\\ \X 2 \ \{ai\ — CL2l)Xi — (022 — a\2)X2 1 , (l.o) 

where A:=detA. For d = 3, the similar case Ilj is defined by ajo = a^o =: a*o, 
i.e., by the two species other than species i having equal growth rates, and by 



Table 1. For d=2, cases when < 2 restrictions yield a Darboux polynomial / and a first integral: \xi\ 1 \x2\ 2 |/| 



case 


restrictions 


DP / 


cofactor of / 


I 


det(ao, a 2 ) = det(ai, ao) = 


none needed 




I'i 


det(aj, aj) = det(ao, a^) = 


none needed 




II 


aio = «20 = 


(an - a 2 i)xi - (a 22 - 012)^2 


an xi + a 2 2 x 2 


H'i 


(X'l j Ob j j — 




Clii X{ 


III 


7-012 = 


aio«20 + a20«n xi + ai a 2 2 x 2 


an xi + a 2 2 x 2 



Table 2. For d = 3, cases when ^ 3 restrictions yield a Darboux polynomial / and a first integral: \x\\ l% \x2\ l2 \x-.i\ l3 \f\ li 



case 


restrictions 




DP / 


cofactor of / 


I 


det(ao, a 2 , 03) = det(ai, ao 


,03) 


none needed 






= det(ai, a 2 , a ) = 


= 






A 


det(aj, aj, a^) = det(ao, a^, 




none needed 






= det(ao, aj,ai) = 


= 






II, 


o-jo = ako =■ a*o, 




( a jj ~ a kj)xj — {akk — ajk)xk 


a*o ~t~ Q« Xi + ajj Xj -\- akk Xk 




®ji — ®>ki — • a*^, 










aiQa^i — a^Qaa 








115 






&i0 ®ii Xi 


an Xi 




Q'jjQ'kk — UjkQ'kj 








Uli 


r 0jk = 0, aji = aki = 




OjOOfeO + O,k0 a jj x j + aj'0«fcfc Xk 


djj Xj + akk Xk 


III' 


7-123 = 0, aio = 020 = 030 = 


: a*o 


{aji aii){aki aii)Xi 


a*o + an xi + a 2 2 x 2 + a 33 x 3 








{®ij ^jj){^ki &ii)Xj 










- {ctik - dkk){aji - au)xk 





CI- 

1' 

o 

a 
«> 
1 

3 

§ 

o" 
a 

o 

c-t- 

r 
1 

1 

to 
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a ji = o-ki =: a *ii i- e -> by their having equal effects on species i. A third condition 
OiO a *i = a*o a ii must also be satisfied. (Note that it will be satisfied if all growth 
rates are zero.) The resulting first integral is 

J I rr . \ a *i( a kj— a jj)( a jk— a kk) \r r .\{ a ii a kk— a ik a ki){ a kj— a jj) I™, I { a ii a jj ~ a ij a ji ) { a jk ~ a kk ) 

A — \ | | | \ tL k I 

x \(a,jj - a kj )xj - (a kk - a jk )x k \ A . (1.4) 

The d = 3 expression (jl.4p reduces to a power of the d = 2 first integral (j 1 . 3|) when 
k) = (3, 1, 2) and the third species decouples: a*j = 0. 
The tables also list restrictions of type III, which are expressed in terms of the 
quantities 

njk '■= {aji — au)(a k j — a,jj)(ai k — a kk ) + (a k i — au)(aij — ajj)(aj k — a kk ). (1.5) 

But the present paper will focus entirely on d = 2 and especially d = 3 systems 
that satisfy the less complicated restrictions of type II, resp. type Ilj. 

It must be stressed that Lotka-Volterra systems with d = 2, 3 and no self- 
interactions (an = for i = 1, . . . , d) do not fit well into the DP framework of tables 
[1] and |2] One example is the classical d = 2 model, which has the two restrictions 
an = 0,22 = and is quite anomalous. It is integrable with first integral 

/=|xi| a20 |x 2 r ai0 [exp(a 2 ixi-ai 22 ; 2 )]=:r[ 3 \fi\ k . (1.6) 

A framework into which this model fits is only now being developed 
( Llibre &: Zhangl 2009). Each factor /, in (j 1 . 6[) satisfies /j = Kifi for some degree-1 



polynomial cofactor Kj and is the refore a DP or generalized DP, sometimes called 
a 'second integral' ( Gorielv 200ll ). To understand (j 1 . 6 1) one needs the concept of 



the multiplicity of an iny ariant line (or algebraic curve) of a polynomial vector field 
(|Christopher et al.ll2007l ). The DP's /i = x u j% 

= X2 are associated to the invariant 



lines x\ =0, rc 2 = 0, which by default are not multiple. The exponential factor f$ 
in (jl.6p comes from the line 'at infinity' which like x\ = and x 2 = is invariant in 
any Lotka-Volterra model but which when an = a 22 = has a nontrivial (double) 
multiplicity. It therefore gives rise to an extra (generalized) DP, namely ^3. 

The extent to which imposing parametric restrictions on a d = 2 Lotka-Volterra 
system can produce single or multiple invariant lines, thereby engendering DP's 



or generalized DP's, is now fully understood (jSchlomiuk &i Vulpd 120101 ). The 
number of invariant lines, counted with multiplicity, can if finite be as large as 
six (including the one at infinity). The line configuration 4.18 of Schlomiuk &: 
Vulpe corresponds to the classical model (an = a 22 =0), and their 4.5 and 4.1 
to the types II and III of table [TJ But an integrability study of their many 
other configurations, most defined by severe parametric restrictions and many 
including multiple invariant lines, remains to be carried out. A d = 3 counterpart 
to their d = 2 analysis, which is lengthy, is not yet available. The many possible 
configurations of invariant planes have not been fully class ified, though partial 



results have been obtained (|Cair6ll2000l : ISaputra et aT]l2010f ) 



There is already a large literature on the integ rability properties of d - 
3 Lotka-Volterra sys t ems without self- i ntera ct ions ( Grammaticos et al. 19901 ; 



Cairo fc LlibTel 120001 : iMoulin Ollagnierl 120011 . 120041 ). They are called ABC 
or A1A2A3 systems, since if an = a 22 = 033 = and aij(ij k a k i^0 for some 
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permutation ijk of 123, permuting species and scaling columns will yield an 
interaction matrix 

/ A 2 1 \ 

A = ( % )? i=1 = 1 A 3 . (1.7) 

V Ax 1 / 



If an ABC system satisfies ^iv^^ = — 1 or Ai = l for some i, and growth rates 
are suitably constrained, a first integral for it can be constructed from degree-1 
DP's. This is because such a system is of type I or type II, as defined in table [2] 
But there are exotic ABC systems with first integrals based on DP's of degree 2 or 
greater, or on generalized DP's. Being tightly restricted parametrically, they are 
not cover ed by table [21 Some of these systems were first obtained by a Painleve 
analysis ( Bountis et al.l 1984]). To some of these exotic ABC systems exotic d = 2 
systems are associated, as any d = 3 system with equal gr owth rates can be reduced 
to a d = 2 system, typically with unequal growth rates ( Cairo et al.ll2003| ). 

A d = 3 system with even tighter parametric res trictions is the May-Leonard 
model of cyclic competition ( Mav &: Leonard] 1 19751 ). This is a system with equal 
growth rates (reducible to zero by a change of variables), and with 



(a, 




'1. 



If, e.g., a + (3 = —1 or a = (3 = 1, then det A = and this system is of type I in 
the sense of table [2] If a = (3 (so that it is fully rather than cyclically symmetric), 
then it is of type II, and specifically it is case Ilj for each i. Hence it is completely 
integrable: for each permutation ijk of 123, the expression (|1.4p will be a first 
integral (cf. IStrelcvn fc Woiciechowski (1988)). Normalized, this is 



I a+1 



ia+1 



-l-2a 



' / |- ; /| ,/ \^k\ \^j %k\ • (l - ^) 

Other first integrals are known dLlibre fc Vallsll201ll : ITudoran fc Oirbanll2012h . 

The relation between the integrability (partial or complete) of a differential 
system and its having the Painleve property (PP), i.e., the property that each 
of its solutions x = x(t) has a one-valued continuation to the complex r-plane 
(without branch points of any order) , is a bit murky. It is widely believed that the 
integration of any system with the PP should reduce to quadratures; or if not, that 
the solutions of the system should be expressible in terms of 'known' functions, 
such as the Painleve transcendents that arise as solutions of second-order ODE's 
with the PP. Like the d = 3 ABC system, the general d = 2 Lotka-Volterra system 
has been subjected to Painleve analyses (IHua et al.lll996l ; iLeach Miritzisll2004l ) . 
It appears that the only such systems with the PP also have DP's from which a 
first integral can be constructed; thus they are at least partially integrable. But a 
Painleve analysis of general d = 3 Lotka-Volterra systems has yet to be performed. 
And until now, no explicit solution involving Painleve transcendents of any d = 3 
Lotka-Volterra system seems to have been published. 



(c) Overview of following results 

In §[2] we illustrate our methods by integrating the above case II of the d = 2 
Lotka-Volterra system with the aid of a 'new time' t. The resulting expressions 
for x\,X2 and r as functions of i, involving an incomplete beta function, are new. 
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This approach is equivalent to that of Gorielv ( 1992h . but unlike him we do not 
posit a fixed form for the new time transformation t = t(r), or make explicit use 
of the first integral (|l,3p . Equal but nonzero growth rates are also treated. 

In §[3] we obtain our central result: the general solution of case Ilj of the d= 3 
Lotka-Volterra system, initially requiring that the growth rates ai0)fl20) a 30 be 
equal. Again the function t = t{r) is not constrained a priori. It turns out to 
satisfy a generalized Schwarzian equation (gSE) , which is integrated with the aid 
of the incomplete beta function. No explicit use is made of the first integral f j 1 . 4 1) . 
Many systems are explicitly solved, including ABC and May-Leonard ones. 

Section [3] re l ates our solution of case-IIj d = 3 systems to the results of 
ICarton-LeBrun ( 19691 ) . who classified all gSE's with the Painleve property. Her 
results yield a classification of equal growth rate case-IIj systems with the 
property: each such has elementary or elliptic general solutions x = x(t). In §|4]we 
express the general solutions of several d = 3 systems with unequal growth rates 
in terms of Painleve transcendents. These may be the first such solutions ever 
obtained. 

In §[5] we summarize the results of §§[2H2 an d make some final remarks. 



2. Two-dimensional integration 

To show how certain Lotka-Volterra systems of the form (jl.ip can be integrated 
parametrically with the aid of a new time variable t, consider the case-II d = 2 
systems of table [1] ones with aio = 020 = 0, i.e., with zero intrinsic growth rates. If 
the interaction matrix A satisfies an 7^ 021 and 022 7^ 0,12, by scaling (redefining) 
the components xi,X2 one can scale the columns of A so that 

IK=( ~ ai 1 ~ a2 ) (2.1) 
^ 1 - a x -a 2 J v ' 

for some a\ , 0,2 ■ The resulting rather symmetric system 

xi = xi \-a\ x\ + (1 - a2)x 2 ] , 

(2.2) 

£2 = x 2 [+(1 - ai)xi - a 2 x 2 ] 

has four invariant lines: the usual x\ = 0, x 2 = and the one at infinity, which 
plays no role here; also x\ — X2 = 0, since if xi = X2 then xi = ±2. Only solutions 
x = x(t) not lying along any of these lines will be considered. By examination, 

I=\xi\ a2 \x 2 \ ai \xi - x 2 \ 1 - ai - a2 (2.3) 

is a first integral: the reciprocal of the standard case-II first integral (|1.3p . 
Define an auxiliary (new time) variable t by 

t:={xi+x 2 )/(x 2 -xi), (2.4) 

so that t = — l,l,oo correspond to the three just-mentioned invariant lines. 
The three t-intervals (— 00, — 1), (—1,1), (l,oo) correspond to sectors in the 
x i, X2-plane lying between the lines. Thus in the first quadrant, the sectors 
< x 2 < xi and < xi < X2 correspond to the i-intervals (—00, —1) and (1, 00). 
To any solution x = x(r) of the system (|2.2p there is associated a function t = t(r), 
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taking values in one of the three t-intervals. The value to taken by t at the time- 
origin (say, at r = 0) determines the t-interval. 

By calculus applied to (|2.4p and the system (|2.2p . t = t(r) satisfies 

*" fl ~ ai + 1 -f^-)=0, (2.5) 



i 2 \t + l t - 1 
which integrates to the hyperlogistic growth law 

i = Ax\t + l\ 1 ~ ai \t-l\ 1 - a2 , (2.6) 
A 7^ being arbitrary. Hence viewed inversely as a function of t, r is given by 

rt 

T = A- 1 x \t' + l\ ai ~ l \t' -l\ a2 ~ l dt'. (2.7) 

-t 

Moreover, it follows by differentiating (|2.4p and exploiting (|2,2p that 

( Xl , X2 ) = ( Xl (t),x 2 (t))=(r^, j^j. (2.8) 

Thus Xi, X2 as well as r can be expressed as functions of the new time t. 

The integral in (|2.7p defines an increasing function of t: an incomplete beta 
function t = B aijO2; f (t) with parameters ai,a2, as is explained in the appendix. 
It maps in an increasing way the t-interval containing to onto some r-interval, 
which may be infinite. The system (|2,2p is thus solved parametrically by 

T = r{t) = A- 1 x B aiia2; t (t), (2.9a) 
(x 1 ,x 2 ) = (x 1 (t),x 2 (t)) = Ax\t + l\ 1 -^\t-l\ 1 ' a2 f^i-p , (2.9b) 

(|2.9b|) coming from (|2.8p . The parameter t is restricted to the relevant t-interval. 

This is a complete integration of the d = 2 system (|2.2p , as the expressions for 
xi,X2 and r as functions of the new time t involve two undetermined constants: 
A and to- (Other than determining the t-interval, the latter merely shifts r.) 
Incomplete beta functions are supported by many software packages, so the 
numerical solution of any case-II d = 2 system is quite easy. Growth rates aio, 020 
that are nonzero but equal can readily be incorporated; see example 12 . 1 1 below. 

One can also eliminate t, and at least formally express x±,x 2 as functions of 
the original time r. Let t = B~ 1 (J2 (r) be any convenient, standardized solution of 

the A = 1 case of the nonlinear ODE (|2.6p . (The function B~ 1 (l2 will map in an 
increasing way some T-interval (r m i n , r max ), which may be infinite, onto one of the 
t-intervals (—00, —1), (— 1, 1), (1, 00); so it will depend on the choice of t-interval.) 
The general solution x = x(r) of (|2.2p is then given by (|2.8p or (|2.9bp with 

t = B- 1 \ a2 {A(r-r )), (2.10) 

which is defined if A(t — tq) lies in (r m i n ,r max ). The solution x = x(t) thus 
contains two free parameters, A and To, and implicitly a choice of t-interval, i.e. 
a choice of sector in the x±, 22-plane. 

For this x = x(r) to be more than formal, an explicit expression for the 
(standardized) inverse incomplete beta function t = B~ 1 1 a2 (r) is needed. As the 
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appendix explains, expressions are available for some a±, 0,2- If the set {1/ai, l/a2, 
1/(1 — a\ —02)} is any of {l,m, — m} (for m a positive integer), {l,oo,oo}, or 
{2,2,oo}, the inverse is expressible in terms of elementary functions; and if it is 
any of {2,4,4}, {2,3,6}, or {3,3,3}, in terms of elliptic ones. In each case the 
inverse extends to a one- valued function on the complex r-plane. Thus x = x(t) 
given by (j2.9b[) and (|2.10p . containing fractional powers, is at most finite- valued. 
Table UH] gives expressions for the (standardized) inverse function t = B~*(r) 

when a = 0, ^, ^, |, 1. Each is one- valued on the r-plane, as stated. It should also 
be noted that if 01,02 are integers which are positive, or satisfy a\CL2 <0 with 
o-\ + 0,2 ^ 0, the inverse will be an algebraic function (and hence, finite-valued 
on the r-plane). In these cases too, x = x(t) will be finite- valued. In any of the 
preceding finite- valued cases, the Lotka-Volterra system (|2.2p has by definition a 
weak form of the Painleve property. This result seems to be new. 

EXAMPLE 2.1. Consider the symmetric system 

±1 = a* xi + xi(-x\ + 3x 2 )/4, 

(2.11) 

X2 = a*o X2 + x 2 (+3 x\ - x 2 )/4, 

where xi,X2>0 and x\ 7^ X2 for simplicity. The state variables x\,X2 are the 
populations of species that grow logistically and also display mutualism: the 
growth of each is made more rapid by the presence of the other. Suppose initially 
that a*o = 0, so that this is the case a\ = <22 = 1/4 of the system (|2,2p . 

The auxiliary variable t = (x\ + X2)/(x 2 — x\) satisfies the interval condition 
t S (1, 00) if (initially and hence subsequently) < x\ < X2; suppose this to be so. 
By table IA~T1 the standardized inverse function for this interval is 

* = B i"/4,i/4( r ) = ^ cq2 + cn ~ 2 K T / 2 )' (2-12) 

where cn is the Jacobian function with parameter m = A; 2 = l/2. The domain 
(r min ,w) is (0,K 1/4 ), where K 1/4 = T(i) 2 /2r(i) « 3.708. Substituting (|2l2| 
into (|2.8p yields the general < x\ < X2 solution, parametrized by A > and tq: 

( Xi ( t ),x 2 (t))=a(-^-, — ) (A(r-r )/2), (2.13) 
' \ 2 cn dn sn cn / ' 

which is defined if A{t — tq) lies in (0, ^1/4). Here such elliptic-function identities 
as sn 2 + cn 2 = 1 and m sn 2 + dn 2 = 1 have been used. For this system the first 
integral (|2.3p is /= |xiX2| 1 ^ 4 |xi — X2\ 1 ^ 2 , and by examination /= \/2A. 

The phase portrait in the < x\ < X2 sector is evident from (|2.13p . 
The endpoints r mm = 0, r max = i^ 1 / 4 , are zeroes of sn, cn respectively. So as 
A(t — To) — > + , t — > 1 + and (x\,X2) diverges to infinity while approaching the 
positive X2-axis (an invariant line). Also as A(t — tq) — > Ky±, t — > 00 and (x%, X2) 

diverges to infinity while approaching the x\ = X2 line (also invariant). These 
divergences are due to the strong mutualism. Although the solution (|2.13p has a 
real r-interval of width A K-^u as its domain, it extends in a one-valued way to 
the complex r-plane. Thus the system (|2.1ip has the Painleve property. 

Incorporating a nonzero common growth rate a*o is straightforward, as noted 
in the introduction. Exponentially modifying the independent and dependent 
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variables yields the general < x\ < solution, parametrized by real T\ and C: 

(sn^ 2dn^\ 
— — - (( e «.o(r-n) _ c)/2 ) (2.14) 

2 cn dn sn cn / ' 

which is defined if e a *°( T-T1 ) — C lies in (0, .K1/4). Compared to the mutualism, 
the role played by any a^o > in causing finite-time blow-up is minor. 

On the complex r-plane the poles of any instance of the solution f|2.13[> lie 
on a square lattice, since sn, cn, dn are doubly periodic; but those of any instance 
of (|2.14p lie on an exponentially stretched lattice. The pole locations are of interest 
because solutions in the complex time domain o f non-integrable syste ms are 
expected to have irregular patterns of singularities ( Bessis &: Chafee 19861 ). 



The just-concluded example was rather special, as setting a\ = 02 = 1/4 leads 
to one-valuedness on the complex r-plane. It must be stressed that for generic 
a%, CL2, owing to the presence of branch points the general complex-domain solution 
(xi(t), X2(t)) of the system (|2.2p will not be one- valued or even finite- valued. The 
distribution of its poles over its many branches remains to be explored. 

However, the example made clear the behaviour of any real-domain solution 
(xi(t),X2(t)) in the forward-time limit; i.e., as r tends to the upper endpoint of 
its interval of definition. In the first quadrant, either the invariant line x\ = X2 
will be an attractor (as in the example), or it will be a repeller and the lines 
x\ = 0, X2 = will be attractors. In population ecology, if a\, 02 > 1 so that the 
two species compete rather than display mutuali sm, these two possibili ties would 
be called competitive coexistence and exclusion ( Mavnard Smith 19741 ). 



The analysis in this section of the case-II d = 2 system assumed for simplicity 
that an 7^ &21 and 022 7^ 012- This can be relaxed. Suppose the growth rates 
are zero (a\o = 020 = 0) and say, that an = a2i =: a*\- Then by examination the 
three usual invariant lines, namely x\ = 0, X2 = and the one at in finity, have 



respective multiplicities 1, 2, 1 in the sense of Christopher et all (12007); and there 



are no others. (This configuration of lines is numbered 4.20 bv lSchlomiuk &: Vulpej 



( 20101 ) .) The line X2 = thus has an extra (generalized) Darboux polynomial 



associated to it. The first integral (j 1 . 3|) is accordingly replaced by 

I = \xi\ a22 I 1 — 0,12 exp(a*ixi/x2), (2.15) 

exhibiting the DP's x±,X2 and the extra generalized DP. It is convenient to choose 
as the new time variable s := x\/x2- The ODE (12. 5p is then replaced by 



(OV2 



s 2 a\2 — 022 v s 



+ 0.1). (2.16) 



This can be integrated to express r = r(s) in terms of an incomplete gamma 
function, or alternatively (if 022 = —012) the error function erf. This sheds light 
on the solution of a particular Lotka-Vo lterra system with an = Q 21 , involving erf 



and erf , that was recently published (jChandrasekar et al.ll2007l . §4.2). 



Integration of three-dimensional Lotka-Volterra systems 



12 



3. Three-dimensional integration 

(a) General case-IIi systems 

It will now be shown that an unexpectedly large family of d = 3 Lotka-Volterra 
systems of the form (jl.ip can be integrated parametrically by the technique of §|2] 
These are case-IIj systems as defined in table |2j specifically, systems with Ojj = 
aki ='■ a*i, so that the species j, k other than i have equal effects on species i. 
Without loss of generality i = 3 and (j, k) = (1,2) will be taken. Initially it will 
be assumed that aio = 020 = a 30 = 0, i.e., that the intrinsic growth rates are zero, 
though including rates that are nonzero but equal is an easy matter. 

Systems in this family include certain ABC and symmetric May-Leonard 
models, which are treated as examples in SS I3l6l and l3lcl Besi des constructing 
explic it solutions, one can determine from powerful results of ICarton-LeBrunl 
(1969) precisely which systems in the family have the Painleve property that each 
solution x = x{t) extends in a one- valued way to the complex r-plane. Generalized 
case-IIj systems which have unequal growth rates will be examined in §|4] 

Suppose «i3 = a23 =: o*3- If the interaction matrix A satisfies 0117^021, 
a 22 7^ &12 and 033 7^ a*3 7^ 0, by scaling (redefining) the components x\, X2, £3 one 
can scale the columns of A so that 





f —a\ 


1 — 0,2 




H 


1 — a\ 


-a 2 






\ h - ai 


62 - d2 


-l/n 



for some 01,02; b%, 62 an d n 7^ 0, — 1. (n = 00 will signify l/n = 0.) The system 
' ±1 = x\ \—a\x\ + (1 - a 2 )x2 + x 3 ] , 

< ±2 = X2 [+(1 - oi)xi - a 2 x 2 + X3] , (3.2) 

k i 3 = x 3 [(61 - a 1 )x 1 + (6 2 - a 2 )x 2 - x 3 /n] , 

which reduces to the d=2 system fj2 . 2[) when 2:3 = 0, has five invariant planes: 
the usual xi =0, X2 = 0, X3 = and the one at infinity, which plays no role here; 
also x\ — X2 = 0, since if x\ = X2 then x\ = xi- Solutions x = x{t) not lying in any 
of these planes are of primary interest. By examination, 

I= \ Xl \^2-b 2 l^pi-ki \ X3 \ \ Xl _ X2 \n(l~a 1 -a 2 )-(l-b 1 -b 2 ) ^ 

is the specialization of the case-113 first integral (|1.4p . Here and below, 

n:={n + l)/n (3.4) 

for convenience; note that with n = 1 meaning n = 00, i.e., l/n = 0. 

Define an auxiliary variable t as in §|2]by 

t:=(xi+x 2 )/(x 2 -x 1 ), (3.5) 

so that t = —1, 1, 00 correspond to the planes x\ = 0, X2 = and x\ — X2 = 0. Thus 
< X2 < x\ and < x\ < X2 correspond to the t-intervals (—00, —1) and (1, 00). 
The function t = t(r) associated to any trajectory x = x(t) takes values in one 
of the intervals (— 00, — 1), (—1,1), (l,oo). On any portion of x = x{t) on which 
t does not change sign, t can be used as an alternative parameter: a new time. 
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By repeatedly differentiating fj3 .5 [) to obtain formulas for t,t,i and i, each 
time exploiting (|3.2p . and eliminating X\,X2,X3 from these formulas, one finds 

2 



t 



+ (n 



+ 



1 - h - 2n(l - ai) 1 - 6 2 - 2n(l - a 2 ) 



t + l 



t-1 



t 



(l-ai)[fei + n(l-ai)] (1- a 2 ) [6 2 + n(l - a 2 )] 



(t + l) 2 (t-1) 2 

(l-ai)[bl+n(l-ai)]+(l-a2)[62+n(l-a2)]+(2-ai-a2)[(l-6i-62)-n(2-ai-ffl2) 

(*+!)(*-!) 



(3.6) 



0, 



which is the d = 3 counterpart of (|2.5p . The nonlinear third-order ODE 
will be called a generalized Schwarzian equation (gSE), since when n= 1/2 (i.e., 
n = —2) its first two terms constitute a Schwarzian derivative, and if additionally 
b\ = ai, 6 2 = a 2 so that its two t/t 2 terms vanish, it becomes a Schwarzian e quation 
of a type familiar from the conformal mapping of triangles (jNeharilll952f ). 

The seemingly complicated gSE (|3.6p . with solution t = t(r), can be integrated 
without difficulty Suppose for simplicity that t > 0. Let t =: i.e., f:=i n , 
and view t, / (rather than r, t) as the independent and dependent variables. By 
applying the chain rule d/dr = i Dt = f l l n Dt where Dt := d/dt, one finds that 



A 2 + 



1 - bi - 2n(l - ai) 1 - 6 2 - 2n(l - a 2 ) 



t + l 



+ 



t - 1 



+ n 



(l-ai)[6i + n(l-ai)] , (1 - a 2 ) [6 2 + n(l - a 2 )] 



+ 



+ 



(3.7) 



(l-ai)[6 1 +n(l-ai)] + (l-a2)[b2+ra(l-a2)]+(2-ai-a2)[(l-bi-&2)-ra(2-a 1 -a2)] 

(t+l)(t-l) 



/ = o 



is the ODE (surprisingly, a linear one) satisfied by f = f(t). This ODE is a so- 
called Papperitz equation, which can be greatly simplified by a substitution of 
a standard type. Substituting / = (t + l)™( 1_ai )(t - l)^ 1 "^) g reduces it to the 
degenerate hypergeometric equation 



1 



DA 9 = 0, 



t+l t-1 

which can be integrated by inspection. In this way one deduces that 



i n = f(t) = \t+l 



in(l— ai) 



t - 1 



\t + l\ n{1 - ai) \t 



1 



in(l— 02) 



in(l— 02) 



i^i + ivT 2 



It' + ll 61 - 1 !^ 



1 



162—1 



(3i 



fit' 



to 



[Ki + K.B^^^^t)], 



(3.9) 



where K\,K% are undetermined constants and to is any convenient point in the 
relevant t-interval. For the incomplete beta function B^ ^.^t), see the appendix. 

Formulas for the inverse function r = r(t) and the alternatively parametrized 
trajectory x = x(t) of the system (|3.2p follow immediately from (|3.9p . The 
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restriction t > will now be dropped: suppose instead that on the portion of 
the trajectory under study, ±i > 0. By integrating dr/dt = i _1 one obtains 



rt 

t = r(i) = t ± 



t ' + l|ai-l|^_l|« a -l|if 1+ X 2 B 6li6a . t0 (t')| 1/(n+1) 1 dt 

to 



/Ml/(n+l)-l ^ 

(3.10a) 



(x 1 ,x 2 ,x 3 ) = (x 1 (t),x 2 (t),x 3 (t)) = (jIj, - (1 - ai)^ - (1 - a 2 )j4rj . 

(3.10b) 

In (|3.10ap . t is restricted to the relevant i-interval (containing any conveniently 
chosen point to) and is further restricted by the requirement that the 
quantity K\ -)- K 2 brt n(t')i within absolute value signs, not change sign. The 
accompanying formulas (13.10bp , which extend the d = 2 formula (|2,8p , are an 
easy exercise. They come by reverting the abovementioned expressions for t, t, t 
as rational functions of xi,x 2 ,x 3 . In practice one would compute t, t in (13.10bp 
by applying i = (dr/dt)" 1 and t = i(d/dt)i to (I3.10al) . 

The parametric solution (r, x) = (r(t),x(t)) displayed in (|3.10p is the central 
result of this paper. It is a complete integration of the d = 3 Lotka-Volterra 
system (|3.2p . as the expressions for r and xi,x 2 ,x$ as functions of t involve 
three undetermined constants: the coefficients Kx,K2, and also tq (which merely 
shifts r). Growth rates aio,a20i°30 that are nonzero but equal can easily be 
incorporated in (|3.10p . as when d = 2. The solution (|3.10p subsumes the d = 2 
solution (|2.9p . This is because if K 2 = then r = r(t) given by (|3.10ap reduces 
to (|2.9ap . and X3 =x^{t) computed from (j3.10b[) will be identically zero. 



EXAMPLE 3.1. Consider a d = 3 Lotka-Volterra system of the form (|3.2p with 
bi = a±, 62 =a>2i so that species 3 is unaffected by species 1 and 2. In this case 
(r, x) = (r(t),x(t)) can be simplified, since if K 2 7^ 0, Eq. (|3.10ap integrates to 

(t\= + {\ K ^ l {\ K ^+ K ^ 1 ,a 2 ;to{t)\ l/{n+1) -iKil 1 /^ 1 )} n finite, 
T TU T ° Xm-H^l + iKz/K^a^M n = oo. 

In terms of a standardized inverse beta function t = B a2 (r) as discussed in the 
appendix (and used in §|2]), this simply says that t = t(r) is of the form 

t - tM -(^(^T + B)^ + C), n finite, 
t ~ t[T) ~\B- 1 \ a2 (eMAr + B) + C), n = oo, [6 ' U) 

for some A^O, B,C. For ai,a 2 for which B~ 1 1 a2 (r) is expressible in terms of 
standard functions, such as in those in table I A, ll this allows x = x(t) when b\ = a%, 
b 2 = a 2 to be similarly expressed, by substituting (|3.12p into (13. 10b|) . 

For many a±,a 2 , such as those in the table, the function t = B~ i 1 fl2 (r) has a 
one- valued extension to the complex r-plane. For such 01,02, if n is an integer 
(n ^ 0, — 1) or n = 00, the system (|3.2p with b\ = a±, b 2 = a 2 will have the Painleve 
property. This is because if n is an integer or 00, t = t(r) in (|3.12p and likewise 
x = x(t) will extend to the r-plane with no branchpoint at r = —B/A, i.e., in a 
one-valued way. But in each such system, the third species decouples. 
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One can also determine which d = 3 Lotka-Volterra systems of the form (|3.2p , 
with no species decoupled, have the Painleve property; and in fact, integrate 
each such explicitly. In the gSE (|3.6p for t = t(r) there are seven terms and 
thus seven coefficients, the first two being 1 and n — 2 = l/n — 1. Suppose the 
remaining five are not determined by a%, 02; 61,62 an d n, as here, but are free 
to vary. For some choices of these five the resulting nonlinear ODE will have 
the P ainleve property. By a mathematically rigorous analysis ICarton-LeBrun 
(1969) determined, classified and tabulated the many possibilities, and integrated 
each resulting ODE. Hence all one needs to do is compare the coefficients in the 
gSE (|3.6p against her tables to determine for which ax, 0,2', 61, 62 and n all solutions 
t = t(r) of the gSE will extend in a one-valued way to the r-plane. For each such 
choice one can compute x = x{t) from (13. 10b[) . and it will extend similarly; thus 
the case-113 system (|3.2p as well as the gSE will have the Painleve property. 

It is not hard to see that for the gSE (|3.6p to have the Painleve property 
the parameter n must be an integer (n 7^ 0,-1 as usual) or 00; else the gSE 
would have complex-domain solutions t = t(r) with movable branchpoints. The 
tables of Carton-LeBrun include infinite families of gSE's in which n is otherwise 
unconstrained, such as ones with 61 = ax, 62 = a<i that were implicitly encountered 
in example 13. II They also include many gSE's with n = 1, 2, 3, 5 or 00. Each can be 
twice integrated to yield a polynomial identity P(t, t; Kx, K 2 ) = with parameters 
Kx,K2, the identities appearing in her table VIII. Any solution t = t(r) is thus 
reduced to quadratures. For nearly all of these gSE's, solutions t = t(r) and hence 
trajectories x = x(t) can be expressed in terms of elliptic functions (lemniscatic 
or equianharmonic) ; for a few, in terms of elementary functions. 

Example 3.2. Consider the d = 3 Lotka-Volterra system 



where for simplicity xx,X2,x 3 > and x\ 7^X2- The state variables Xx,x 2 ,x 3 are 
the populations of three species with common intrinsic growth rate a*o- If a*o > 
then in isolation species 3, though not 1 or 2, grows logistically. Species 1 and 2 
display mutualism, as do 1 and 3; and species 2 preys on species 3. 

Since a nonzero a*o is easy to incorporate, set a*o = 0. System (|3.13p is then 
the case (a 1 ,a 2 ; 61, 6 2 ; n) = (0, 0; 1, -2; 1) of (pT2]t . With these values the gSE (1331) 
for t = t(r) is a variant of Carton-LeBrun's n = 1 case II. Integrating twice yields 



where Kx, K2 are undetermined constants. (As i n =i 2 = f(t), the right side is the 
general solution of the Papperitz equation (|3.7p .) Integrating once more yields 



' ±1 = a* xx + x\{x2 + x 3 ), 
< x 2 = a* x 2 + x 2 (xx + x 3 ), 
, X3 = a* x 3 + x 3 (xx -2x 2 - x 3 ) 



(3.13) 



f = Kx(t 2 -l) 2 + K 2 {t + l) 2 



(3.14) 




(3.15) 



as the general solution of the gSE. Here c\, C2 are constrained by c\ 
t can be replaced by At + B, for any A 7^ and B. 




1 and 
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Substituting ff37T5]) into fj3.10b[) yields 

(x 1 ,X 2 ,X 3 ,)(r) = 

( C2 sinr 2c2 sinr — C2 + (1 + c\) cos r 

\ 1 — (ci — C2 cos t) ' 1 — (ci — C2 cos r) 2 ' [1 + (ci — C2 cos r)] sin T 

as the general solution x = x{t) of the system (|3.13j) with 0*0=0, it being 
understood that r can be replaced by At + B, with x scaled by A. There are 
three degrees of freedom: A, B and the point on the hyperbola c\ — c| = 1. 

The first integral I of fj3.3[) specializes to x\\x^ x 3 \, and by examination 1 = 
A 2 x |ci + l|/2. The existence of a explicit general solution facilitates the study 
of global dynamics. 

EXAMPLE 3.3. Consider the d = 3 Lotka-Volterra system 

' XX = X\ \-\ X X + (1 + })X 2 + X 3 ] , 
< X2 = X 2 [+(1-^)XX + ^X2+X 3 ], (3.17) 
k X 3 = X 3 X X - % X 2 - X 3 ] , 

where r ^ 0, and for simplicity x\, x 2 , x 3 > and x\ ^ x 2 - The system fj3. 13[) is the 
case (oi,a 2 ; 6i, 6 2 ; n) = (£, J, ; 1) of (|5^ . With these values the gSE (I3T61) 
for i = t(r) is a variant of Carton-LeBrun's case e.IV. The gSE can be integrated 
twice to yield t = (u r + l)/(u r — 1), where u is a solution of 

u 2 =K lU A -K 2 (3.18) 

and Kx,K2 are undetermined constants. (To confirm this, it is easiest not 
to integrate the gSE itself, but rather the Papperitz equation (|3.7p or 
preferably the degenerate hypergeometric equation (|3.8p . which are equivalent.) 
By integrating (|3.18p once more, one deduces that u oc cn(r) and thus that 

, xi, v C en ~\~ 1 . . , . 

t = {xx + x 2 )/(x2-x l )= A t). (3.19) 

o cn — i 

Here cn is the Temniscatic' Jacobian elliptic function with modular parameter 
m = k 2 = 1/2, and C is arbitrary; and as above, r can be replaced by At + B. 
Substituting the expression (|3.19p into (I3.10bp yields 

w . / rsndn Crsndncn 1 " -1 cn 3 \ . . 
( Xi ,x 2 ,x 3 )(t)=( , - , — - r 3.20 

\(Gcn r — ljcn C cn r — 1 sndny 

as the general solution x = x{t) of the system (|3.17p . it being understood that 
r can be replaced by At + B, with x scaled by A. The first integral I of fj3. 3[) 
specializes to \xi/x2\ 2 ^ r \x 3 (xi — X2)\, and by examination I = A 2 x |r| |C| _2//r . If 
r is an integer then (|3.20p like (|3.16p will be one-valued on the complex r-plane; 
thus the system (|3.17p like (|3.13p will have the Painleve property. 
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(b) ABC systems 

An interesting application of the results of the last subsection is to a well-known 
family of ABC systems, which as explained in S fllfel are (generic) d = 3 Lotka- 
Volterra systems without self-interactions. Any ABC system, with interaction 
matrix A= (oij)ij=i gi ven by (|1.7p . is of the form 

yi = yi(A 2 2/2 + 2/3), 

m = 2/2(^3 2/3 + vx), (3.21) 
m = yz{Myi + 2/2)- 

It is known that when one of Ai,A 2 , A3 equals unity the ABC system is completely 
integrable in that a pair of first integra l s exis ts, one of them being of the Darboux 
polynomial type ( Grammaticos et al.l Il990l ). But no general solution of ABC 
systems of this type has previously been published. From our point of view the 
integr ability comes from an ABC system with A{ = 1 being a case-IIj system as 
defined in table [2] This facilitates a complete integration. 

Suppose without loss of generality that A3 = 1, and also that A 2 7^ 0. In terms 
of a scaled state vector (x%, x 2 , X3) := (yi, A 2 y 2 ,y 3 ) the system becomes 

■x 1 (x 2 +x 3 ), 

-x 2 (xi +23), 



Xl 



x 2 



X3 



(3.22) 



X2,(Ai xi + A 2 l x 2 ) 



which is the (aj., a 2 ; bi,b 2 ;n) = (0, 0; Ai, A% ; 00) case of the case-Ii3 system (|3.2p . 
An integration of (|3.22p therefore follows at once from the complete 
integration (|3.10p . (As above, solutions that do not lie in any of the invariant 
planes x\ = 0, x 2 = 0, x\ — x 2 = are of primary interest.) 

The new time t := {xi + x 2 )/(x 2 — x±) is used as the parameter. Along any 
segment of a trajectory x = x(t) on which ±(t = dt/dr) > 0, the variables r and 
x = (xi, x 2 , x$) are expressed as functions of t by 



r(t) = t ± 



to 



[(f) 2 -1] [Ki + K 2 B 



(xi(t),x 2 (t),x 3 (t)) 



t + r t-r t t + i t-i 



(3.23a) 



(3.23b) 



Here tq is the initial time, at which t = tQ. The parameter t is restricted to 
whichever of (— 00, — 1), ( — 1,1), (l,oo) contains to. It may be further restricted 
by the additional condition that K\ + K 2 B^ A' 1 -^^')' within the absolute value 

signs, not change sign over the integration interval. The derivatives t, t in (j3.23b[) 
are computed by applying t = (dr/di) -1 and t = i(d/dt)i to (|3.23ap . 

The parametric general solution (|3.23p is a complete integration of the ABC 
system with A3 = 1, since the expressions for r and xi,x 2 ,x% involve three 
undetermined constants: Ki,K 2 and also To (which merely shifts r). For many 
choices of Ai,A 2 the expressions for r and xi,x 2 ,x 3 in terms of t can be 
expressed using elementary functions. This is because for many choices of b\, b 2 the 
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incomplete beta function ^b 1 ,b 2 ;t (t) in (|3.23ap can be so expressed. The evaluation 
of t and x\,x 2 ,x% is especially easy if ^.^(i) is a rational function of t. 

Example 3.4. Consider the case (A\, A 2 , A s ) = (1, 1, 1) of the ABC 
system (|3.2ip and hence of (|3.22p . for which (ai, a 2 ; bi, b 2 ; n) = (0, 0; 1, 1; 00). 
When t > 1, as is the case for any trajectory x = x{t) in the sector < x\ < x 2 , 
the definition (jA.2|) implies that 'Bi t i-t (t) equals t plus a constant. Substituting 
into (|3.23p and evaluating the integral and derivatives yields 

r(t) = "c^r + w^c) + WTcy (3 - 24a) 

(xi, x 2 , x 3 )(t) = ((t - l)(t + C), (i + l)(t + C), t 2 - 1) , (3.24b) 

where C is a constant; and one can also replace r by At + £? and scale (x\, x 2 , x$) 
by A. Thus there are three undetermined constants in all. 

The solution x = x{t) in the sector is generated parametrically by allowing t to 
range over (1, 00), or perhaps (depending on C) a proper sub-interval. (The latter 
is because t > 1 covers > x\ > x 2 as well as < x\ < x 2 .) The first integral I 
of (|3.3p specializes to |x^~ 1 x 2 ~ 1 X3(xi — a?2) 2 | 5 an d by examination I = 4|t4|. 

Trajectories in other sectors, such as <x 2 < x\, are generated similarly. The 
general solution f|3.24|) resembles but is more ex plicit than a solution of the ABC 
system with (A x , A 2 , A 3 ) = (1, 1, 1) obtained bv kbenda et all (]200ll . § 7). 



For any ABC system of the form (|3.22p . the new time t = t(r) satisfies the 
generalized Schwarzian equation (|3.6p with (qj , q 2 ; fri , b 2 ; n) = (0, 0; A-\ , A^ 1 ; 00) . 
(That n = 00 implies n = 1.) Several of the gSE's listed bv ICarton-LeBrunn 19691 ) 



as having the Painleve property, in that all solutions t = t(r) extend to a one- 
valued way to the complex plane, have n = 00 and a± = a 2 = 0. They therefore 
give rise to ABC systems with ^3 = 1 that have the Painleve property. 

EXAMPLE 3.5. Consider the case (A±, A 2 , A3) = (-1/2, -2, 1) of the 
ABC system (|3.2ip and hence of (|3.22p . for which (01, a 2 ; bx, b 2 ; n) = 
(0,0; -1/2, -1/2; 00). With these values the gSE (l3T6|) satisfied by t = t(r) is a 
variant of Carton-LeBrun's case e.IV.3. The solution of this gSE is 

t( T ) = [{C tanh r) + (C tanh r)" 1 ] / 2, (3.25) 

where C 7^ is arbitrary and r can be replaced by At + B. (This is best shown 
by integrating the hypergeometric equation (|3.8p . though one could also integrate 
Eq. (|3.23ap . as in the last example.) Substituting f|3. 25[) into (|3.23b[) gives 

(xi,x 2 ,x 3 )(t) = (3.26) 
— cosh + C sinh — cosh — C sinh 2(1 — C 2 ) sinhcosh\ 

1 ( 7" ) 

sinh cosh (cosh + C sinh) ' sinh cosh (cosh — C sinh) ' C 2 sinh 2 — cosh 2 / 

where again r can be replaced by At + B, in which case (xi,x 2 ,xz) is scaled 
by A. Thus there are three undetermined constants in all. 

Equation (|3.26p provides a complete integration of the ABC system with 
(Ax, A 2 ,A 3 ) = (-1/2, -2, 1). If -1< C < 1 the trajectory defined by (HI2TD with 
r < lies in the positive orthant; and if C = it lies in the plane x\ — x 2 = 0. 
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If C = ±1 it lies in the plane x 3 = 0. The first integral / of (|3.3p specializes to 
\x\X2\ l ^ 2 \x^/{xi — £2)1, and by examination /= \A\ x |(C 2 — 1)/2C\. 

This system is of much interest, as it is the A 3 = l member of a family of ABC 
systems (with A x = -1/(43 + 1) and 42 = -(^3 + l)/^3 for 4s 7^ 0, -1), each o f 
which has the Painleve property and is completely integ rable (Bountis et al.lll984l). 
Each has a first integral that is quadratic in xi, x%, x 3 ( Strelcvn fc Woiciechow ski 



Il988h . If A3 = 1 this extra first integral is J = (x\ — 22) + ^ x ?,{ x i + x 2 + X3), and 
by examination J = A 2 x 16, irrespective of the integration constant C. 

Bountis et al. remark that the member with (Ax, A2, A3) = (—1/2, —2, 1) can 
be integrated with the aid of elliptic functions; but as Eq. (|3.26p reveals, sinh and 
cosh suffice. This explicit general solution also makes manifest the one-valuedness 
of the continuation of x = x(t) to the complex r-plane. 

(c) Symmetric May-Leonard systems 



A particularly striking application of the results of § l3lal is to the integration 
of certain May-Leonard systems. As was explained in § \l\b\ the three species 
in any May-Leonard system compete cyclically, with a, f3 being the strength 
of clockwise and counter-clockwise interactions. The systems of interest here 
are those with a = J3, for which the species are equivalent; they are completely 
integrable. An explicit general solution will be obtained for the first time: see Eqs. 
(|3.29p - (|3.30p below. The system with a = (3 = — 1, which is mutualistic rather 
than competitive, will be shown to have the Painleve proper ty. May-Leonard 
systems have been analysed fro m the Painleve point of view (IBountis fc Segur 
ll982l:ISachdev fc Ramananlll997l) and multiple first integrals have been discovered 
(|Llibre fc Valid l201ll : iTudoran fc Girbanll2012f ). But the results below are much 
more explicit. 

Suppose that a\o = 020 = 030 = 0, i.e., that the intrinsic growth rates are zero, 
since including rates that are nonzero but equal is easy. The system (jl.ip . with 
the May-Leonard interaction matrix (j 1 . 8 1) and a = /3, is then of the form 



Vi = Ui{-yi -oiy 2 -ay 3 ) 



V2 



yzi-ocyi 



y2 ~ay 3 ) 

o-yi- yz) 



(3.27) 



Suppose that a 7^ 0,1. Then in terms of a scaled state vector {x\, X2, X3) := 
((1 — a)yi, (1 — a)y2, —ay^), the system becomes 



Xl : 

X2 



xii-j^xi -j^x 2 + x 3 ) 
x 2 (- a - 1 



X3 =X 3 (- 



1-a 



l-a 



Xl 



l-a 



X2 + X 3 ) 



(3.28) 



xi - tz^ x 2 + 



X3), 
1. 1 v 

. .1.1. 



which is the (ai, 02; 61, ^2; n ) = ((1 — o)^ 1 , (1— a) -1 ; 1,1; —a) case of the 
case-113 system (|3.2p . An integration scheme therefore follows at once from the 
complete integration (|3.10p . As above, it will yield all solutions that do not lie in 
any of the invariant planes xi =0, X2 = 0, xi — X2 = 0. 
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The new time variable t := (x\ + x 2 )/{x 2 — %i) plays its familiar role. Along 
any segment of a trajectory x = x{t) on which ±(i = di/dr) > 0, the function 
t = t{r) satisfies 

i = ±\(t 2 - l){K 1 t + K 2 )\ a/(a ' 1) (3.29) 

where Ki,K 2 are constants. This is a specialization of (|3.10ap . since according to 
the definition (jA.2|) . Bi j i ; t (t) equals t plus a constant. After obtaining t = t(r) 
by integrating Eq. (|3.29p . one would substitute it into (13.10bp . i.e., into 

ta(r),*,W.„ W ) - ±- r j' + T ^ (±- + ^)) , (3.30) 

to obtain the components of x = x(r). This scheme is readily carried out, either 
numerica lly or (for certai n a) symbolically. Thus one can, inter alia, confirm the 
results of Ble et all ( 2013f ) on the global dynamics of a = [3 May-Leonard systems. 



EXAMPLE 3.6. Suppose a = m/(m + 1), for m a positive integer. Then the 
exponent a/(a— 1) in (|3.29p equals — m, and dr/dt = i~ 1 is a degree-3m 
polynomial in t. Hence r is a degree-(3m + 1) polynomial, and its inverse t = t(r) 
is algebraic. By substituting this function into (|3.30p one deduces that on the 
complex t- plane, each of xi,x 2 ,X3 is also an algebraic (finite-branched) function. 

Hence when a = f3 = m/(m + 1), the May-Leonard system has a weak form of 
the Painleve property. This result appears to be new. The simplest case is when 
a = 1/2 (corresponding tom=l), when r is quartic in t. By applying the quartic 
formula one can express t = t(r) and hence each X{ = Xj(r) in terms of radicals. 

Example 3.7. Let a = /3 = — 1, so that (oi, 02561, b 2 ; n) = (1/2, 1/2; 1, 1; 1). 
When a = — 1 the solution of the ODE (|3.29p can be expressed using the 
Weierstrassian elliptic function p(r) = p(r; g 2 , 53), with p 2 = 4p 3 — g 2 p — 93- One 
notes that the gSE (|3.6p . also satisfied by t = t(r), is an ODE listed by 
Carton-LeBrunl (1969), as her n = l case XIII. 4. Regardless of which ODE one 



integrates, one finds (with C a constant of integration) the general solution 

t(r) = p(r; 52,53) + C, (3.31a) 

92 = 4(3C 2 + 1), g 3 = 8C(C 2 -l), (3.31b) 

where r can be replaced by At + B. By substituting into (|3.30p and using the 
identity p = 6p 2 — g 2 /2 one obtains the general solution x = x(t): 

( V ^ ( P P 2[(p + C) 2 -l] ^ 

(xi,x 2 ,x 3 )(r) = , 7—7^ — -, — -. L , (3.32) 

yp+C+lp+C-1 p y 

where r can be replaced by At + B, with x scaled by A. This being meromorphic 
on the r-plane, the a = /3 = — 1 May-Leonard system has the Painleve property. 

Tha t this system can be integrated using elliptic functions was noted by 
Brenid (jl988l ). but the explicit solution (1332]) is new. It is easily checked that 



the functions xi(x 2 — 2x3), x 2 (xi — 2x3) and x 3 (xi — x 2 ) are first integrals: they 
are time- independent, depending only on C and A. They are essentially the first 
integrals introduced as Ii,I 2 ,I 3 in (|1.9p above (where x is to be read as y) . S t arting 
fro m (13.321). one can confirm th at the first integrals found bv lLlibre Vallsl ( 2011 ) 
and lTudoran fc Girbanl (|2012l ) are also time-independent. 
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For the general solution x = x{t) of a version of the a = (3 = — 1 May-Leonard 
system with a single nonzero growth rate, see example l4,ll below. 



4. Unequal growth rates and Painleve transcendents 

Section [3] treated d = 3 Lotka-Volterra systems classified as 'case IL' according to 
table 12 with zero growth rates. (The extension to nonzero but equal growth rates 
is easy as was noted.) Several d = 3 systems with unequal growth rates, which have 
the Painleve property, will now be completely integrated. They are deformations 
of systems treated in §[3] but do not fit into the framework of table [2j This is 
because their first integrals, based on Darboux polynomials, are not conserved: 
they depend exponentially on r. Their solutions x = x(t) turn out to involve 
Painleve transcendents, which for Lotka-Volterra systems is a novelty. 
Consider the system 

' x\ = x\ \-a\X\ + (1 - a 2 )x 2 + x 3 ] , 

x 2 = X2 [+(1 - ai)xi - a 2 x 2 + x 3 ] , (4.1) 

. x 3 = x 3 [A + (bi - a x )xi + (b 2 - a 2 )x 2 - x 3 ] , 



which is the n = 1 specialization of the case-113 system (|3.2p . modified to include 
a growth rate A for species 3. Up to scaling, this is the generic system in which 
species 1,2 have equal effects on species 3, species 3 has an equal but negative 
effect on itself, and species 1,2 have zero growth rates. As in §§[2]and[3l define an 
auxiliary variable t by t = (xi + x 2 )/{x 2 — x\). By calculus one deduces that 

(x1.x2.x3) (r) = (J!, jij, f - (1 - oij^j - (1 - a 2 )j^) (r), (4.2) 

which previously appeared as Eq. (]3.10bp . 

The new time t = t(r) satisfies an extended version of the generalized 

Schwarzian equation (|3.6p . If the left side of (|3.6p is abbreviated as gS(ai,a2; 
bi,b 2 ;n), by tedious elimination one deduces that 



gS(ai,a 2 ;6i,6 2 ;l) + A 



1 — a± 1 — a 2 \ 1 1 



t + 1 t - 1 J t i 3 



0. (4.3) 



It is straightforward to rewrite this nonlinear third-order ODE in the form 
J — A J = 0, where 

J ■= V + L (A A) 

(t + l)l-2ai+6i(i _ l)l-2a 2 +6 2 V ' 

can be viewed as a time-dependent first integral. The existence of such a 
representation is unsurprising: by exploiting (|4.2p . one can see that J is the 
familiar DP-based first integral f j 3 . 3 [) of §|3j written in terms of t, t, t (and with 
absolute value signs omitted, so that in general, branch choices must be made). 

To deal with the exponential dependence on r, introduce a new independent 
variable z:=e Xl ~/ c for some c^O, and denote d/dz by a prime. The statement 
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J = Joe can be written as a nonlinear second-order ODE for t = t(z), i.e. 

, = (Izsl + (0 » - £ + g £± g^jb D-^ , (4 . 5) 

where K = c 2 Jq/X 2 is a constant of integration. Equation (|4,5p is a 'proto- 
Painleve' OD E, in tha t for certain a±, 61, b%; c, it defines a Painleve 
transcendent (lined 11927 ) . This is obscured by its singular points being t = 
— l,oo, +1, while the defining ODE's for the transcendents (traditional; due to 
Painleve and Gambier) have singular points w = 0, 1, 00. An equivalent ODE is 

V w> to - 1 / 2 z 2 ~ c 

where tu = (t + l)/(i — 1) = x%/xi is a new dependent variable, with t = (w + 1)/ 
(u; — 1); and by definition 03 := 1 — ai — 02 and 63 := 1 — b\ — 62- 

Example 4.1. For certain (ai, 02; 61, 62; c), the ODE (j4.6|) for to = is or is 
reducible to a Painleve-III or Painleve- V equation. The Pni and Py equations have 
parameters a, (3,^,5; and for N = III or V, w^(z) = w^(a, /?, 7, 5; z) will denote 
any solution of the respective equation. The possibilities include 

(01,03; 61,625 c) = (0,0; ^-| ; 2) : w = w m (0, 0, K\ 0; z) r , 

= (0, 0; -§, -i; |): w = V n (if', 0, 0, -if'; z) 3 , 

= (0, 0; -|, -|; 1) : u; = iwm^', A", 0, 0; z) 2 , 

= (0, i; -|, 1; 1) : w = i?Kn(if' , A', 0, 0; z)% 

= (±,±;1,1;1) : w = wy(0, 0, K' , 0; z), 

where r 7^ is arbitrary, R[w] := —&w/(w — l) 2 ; and A' oc if is free in each case. 
These are Carton-LeBrun's n = l cases e.IV, XI.l, XIII. 1, XIII. 3, XIII.4. By 
substituting t=[w + X)/{w — 1) into (|4.2p . one obtains (xi, £2, £3) in terms of 
the transcendent w^{z = e Xr ^ c ) and its derivatives. 

It should be noted that the last of these complete integrations, in terms of wy, 
is of a deformed (A 7^ 0) version of the a = f3 = —1 May-Leonard system, which 
was solved (with A = 0) in example 13.71 

The preceding was stimulated by later ( 1970f ) results of Carton-LeBrunl . who 
studied and integrated many gSE's of an extended type that includes Eq. (|4.3p . 



5. Summary and final remarks 

We have shown [in Eq. (|3.10p ] how to construct the general solution x = x(t) 
of any d = 3 Lotka-Volterra system in which species j, k have equal effects on 
species i, and the three species have equal growth rates. Such systems had not 
previously been integrated, despite extensive work. The constructed solution is 
parametric, with r, x expressed with the aid of the incomplete beta function as 
functions of the new time variable t := (xj + Xk)/{xk ~ x j)- If the system has the 
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Painleve property, the new time is typically not needed: as a function of r, the 
system state x can be expressed in terms of elementary or elliptic functions. 

From the complete integration of any Lotka-Volterra system of this type, one 
can derive first integrals. One is the DP-based first integral of (|1.4p . the constancy 
of which defines type Ilj. A n additional firs t integral involving a qua dratu r e was 
found in the ABC case by iGorielvl ([19921 ). and more generally by I Gaol (j2000, 
theorem 5). But they did not construct general solutions. It seems difficult to go 
from a pair of first integrals to a general solution, rather than the reverse. 

The fully symmetric (a = 0) May-Leonard model was integrated, and the a = 
j3 = — 1 case was treated in examples 13 . 71 and 14 . 1 1 If the growth rates are equal, x is 
expressed in terms of the Weierstrassian elliptic function p(r); and if a single rate 
(A) is nonzero, the Painleve transcen dent roy(e^). This case, w ith equal growth 
rates, was studied geometrically by iTudoran &; Girbanl (|2012l ). but the elliptic 
general solution is new. The appe a rance of a P ainleve function is also novel. 

The results of Carton-LeBrunl ( 19691 . Il970l ) on certain nonlinear third-order 
ODE's with the Painleve property pr oved inva l uable. They are not well known, 
though they have been re- worked by ICosgrovel (jl997l . Appendix C). In many of 
our examples they facilitated the integration of the ODE satisfied by the new time 
t = t(r). With further effort, they could probably be made to yield a classification 
of all type-II, Lotka-Volterra systems with the Painleve property. 

Our techniques may be useful in constructing solutions of other small-<i 
quadratic dynamical systems, of the sort reviewed in the introduction. A related 
inverse problem is also of interest. Many functions of time r have Lotka-Volterra 
representations, i.e., can be generate d from solutions x = x(t) of systems of Lotka- 
Volterra type. (This is stressed by iPeschel fc Mende ( 19861 ). who even mention 
hardware implementations.) Painleve functions are included, as §H] made clear. 
Along this line, the quadratic represent ability of solutions of further higher-order 
ODE's of Painleve type, such as Chazy equations, will be explored elsewhere. 



Appendix: Incomplete beta functions and their inverses 



In a normalization used here, the incomplete beta function 3 ai a2 -t (t), for (real) 
indices 01,02, (real) endpoints —1,1 and a (real) basepoint to ^ —1,1, is given by 



li' + l 



|ii — 1 



t'-l 



\a,2— 1 



dt'. 



(A.l) 



to 



The basepoint to hes in one of the intervals (— oo, — 1), (—1,1), (l,oo), and by 
convention r = B aijO2 .£ (t) is defined only for t in that interval. It is an increasing 
function of t. If a\, a>2 are positive integers it is a polynomial function; and if a±, 02 
are integers with a\a2 < and a\ + a<i ^ 0, it is a rational function. In general it 
is a 'special' (higher transcendental) function. If a\,a<i >0 it can be expressed 
in terms of the traditionally normalized incomplete beta function, which is 



B i (ai,a 2 ) 



t ai -\l-t) a2 - l dt 



(A.2) 
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for t G (0, 1). Specifically, 

{ B (m)/(t-i) (ai, 1 - oi - o 2 ) + C, t G (-oo, -1), 
B (t+ i)/ 2 (ax, a 2 ) + C, t G (-1,1), 

B (t-i)/(m) («2, 1 - ai - a 2 ) + C, t G (1, oo), 

(A.3) 

with C = C ai O2 -t chosen so that B aija2; t (io) = 0. The traditional function 
B£(oi,a2) is supported by many software p ackages, and many identities, 



expansions and approximations for it are known (lDutkalll981[ ). In particular 



B i (a 1 ,a 2 ) = a^ 1 i^(l-i) a2 2 F 1 {a 1 + a 2 ,l; a x + 1; t), (A.4) 

where 2 F\ is the Gauss hyperge ometric function. Thus for many choices of a\,a 2 , 
closed-form expressions exist ( Prudnikov et all Il990. §7. 3), A few have been 



derived in the context of hyperlogistic population growth (iBlumberg 1968). 

The inverse incomplete beta function t = B~ 1 1 a2 . 4o (r) is also an increasing 
function of its argument, and is algebraic if a±, a 2 are positive integers, or integers 
with a\a 2 < and a\ + a 2 ^ 0. This inverse function is defined on some interval 
containing r = 0, and satisfies (with t = dt/dr) the hyperlogistic growth equation 

i=\t + l\ 1 - ai \t-l\ 1 - a2 , (A.5) 

with the initial condition t(r = 0) = B~^ art (0) = to- It maps the r-interval onto 
whichever of (— oo, —1), (— 1, 1), (1, oo) contains to. It is useful to write 

K!,ar,t ^)=^a 1 \a 2 (r-ro), (A.6) 

where t = B~ 1 a2 (r) is any convenient, standardized solution of (1A.5|) that maps 
some r-interval (r m i n ,r max ) onto the t-interval containing to- The time-origin ro 
is determined by the condition that B" 1 ( — ro) =to- 

The function r = B aij(l2; t (t) can be continued analytically from its real interval 
of definition to the upper half of the complex t-plane. (The continuation is given 
by (lA.ip without absolute value signs, multiplied by an overa ll phase factor.) 
The continuation is a Schwarzian triangle function ( Neharilll952 ) that performs a 



Schwarz-Christoffel transformation: provided oi, a 2 , 1 — at — a 2 are non-negative, 
it conformally maps the upper half of the complex t-plane to the interior of some 
triangle AABC in the r-plane with vertex angles Tr(a\, a 2 , 1 — a\ — a 2 ), the points 
— 1, l,oo on the real t-axis being taken to the vertices A,B,C. Many conformal 
mapping functions of this type, which are essentially incomp l ete b eta functions, 
can be found in the catalogue of von Koppenfels &; Stallmannl ( 1959f ). In each case 



the inverse function takes the interior of AABC to the upper-half t-plane. 

This inverse t = B~ 1 (l2 . to (r) can sometimes be given in closed form; e.g., when 
the unordered set {1/ai, l/a 2 , 1/(1 — a\ — a 2 )} is any of {2,4,4}, {2,3,6} or 
{3, 3, 3}. In these cases the inverse can be continued to the entire complex r-plane 
from its real interval of definition, and furthermore from AABC, by applying 
the Schwarz reflection principle: reflecting repeatedly through the sides of the 
triangle. The resulting one-valued function of r is elliptic (i.e. doubly periodic) 
and can be given explicitly. The cases {l,m, —m} (for m any positive integer), 
{l,oo,oo} and {2,2,oo} are similar, but yield inverse functions t = B~ 1 a2 . io (r) 
that are elementary rather than elliptic. 
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Table A.l. Five cases when t = t(r) = B ai 1 , Q2 (r) is expressible in closed form and can be extended 





to a one- valued function on 


]he complex r-plane. 




a 


t = t(r) = B- L a (r), 


T £ ("7~mm; ''"max) 






to G (— 1, 1) subcase 


±t a £ (l,oo) 


subcase 



i 

4 
1 

I 
2 
1 


tanhr, rG (—00,00) 
V2[sndn](T/V2), r e (-K 1/4 , K 1/4 ) 
-81p/[2(9p + l) 2 ], re(-K 1/3 ,K 1/3 ) 
sinr, r G (-7r/2, 7r/2) 

re (-1,1) 


— coth r, 

±i[cn 2 + cn- 2 ](r/2), 
±1-8/(270 ±2), 
± cosh r, 

±1 + T, 


±re(-oo,0) 
±re(0,K 1/4 ) 

±T 6(0,^1/3) 
±T G (0,oo) 
±r G (0, 00) 



Thus there are many 'conformally distinguished' choices for the unordered set 
{1/ai, 1/02}) each of which yields an explicit formula for the inverse function t = 
B~ i 1 a2 . <o (r). The elementary cases are {l,m}, {1,—m}, {m,—m}, {l,oo}, {2,2}, 
{2,oo}, {00,00}, and the elliptic ones are {2,3}, {2,4}, {2,6}, {3,3}, {3,6}, 
{4,4}. In each of these cases any standardized inverse function i = B~ 1 a2 (r) of 
the type used in fjA.6[) . which has a real interval of definition (r m i n ,r max ), can 
optionally be continued to a one-valued function on the r-plane. 

Table IA.1I gives such a standardized function t = t(r) = B~„(t) for a = 

0, |, |, 1. These are the conformally distinguished cases with a\ = a 2 = : a, and 
are of particular interest (the ones with aq 7^02 are left to the reader). In each 
case there are two subcases: to G (— 1, 1) and ±to £ (l,oo). The function t = t(r) 
maps r 6 (r m i n , Tmax) onto t G (—1, 1), resp. ±t G (1, 00), and by (jA.5|) it satisfies 

^(l-t 2 ) 1 "", resp. i=(f-l) 1 - a . (A.7) 

In the table each t = t{r) has been chosen to satisfy t(r = 0) = 0, resp. t{r = 0) = 
±1 (except when a = 0). Each is increasing on its interval of definition (T m j n , T max ). 

The elementary functions appearing in the table (for a = 0,^,1) follow by 
inspection: by integrating the ODE's (jA77j) . The elliptic solutions of these ODE's 
(for a = ^)|) are less obvious, but can readily be derived with the aid of 
elliptic function identities from func tions that are known to map certain tria ngles 
conformally to the upper half-plane ()Koberlll957l :[S ansone fc Gerrets en 1969|). For 
a = \ the Jacobian elliptic functions sn, cn, dn that appear are 'lemniscatic': their 
common modular parameter m = k 2 equals ^ . For a = | the Weierstrassian elliptic 
function p(r) = p(r; g2, ^3) that appears is 'equianharmonic': its parameter 52 
equals zero. Also, its parameter 53 equals — 4/3 6 = —4/729. It satisfies the usual 
Weierstrassian ODE p 2 = 4p 3 — g 2 p — 53, i.e., p 2 = 4(p 3 + 3~ 6 ). 

For each function t = t(r) = B~^(r) appearing in the table, the interval of 
definition r G (r m i n , r max ) is easily computed from (lA.ip . It can be expressed 
in terms of the traditionally normalized complete beta function B(ai,ci2):= 
r(ai)r(a2)/r(ai + 02). In the subcases to £ ( — 1> 1) an d i^o £ (1) °°); it is 

Te(-K a ,K a ), resp. ± t G (0, K a ), (A.8) 
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where 



ir(a)r(i) r , . . ir(o)r(i 



:= | B(a, |) = - Y^prp res P- ^a:=^B(a, 2 ^, - 2 ^ 

It is understood that K a is infinite if a ^ 0. For instance, Kq,Kij2,K\ are 
oo, 7r/2, 1, as shown. Also, in the ±io 6 (1, oo) case the interval ±r 6 (0, X a ) must 
be replaced by ±r £ (— oo, 0) if a ^ 0, and interpreted as ±r E (0, oo) if a ^ ^. 

The values of a in the table are not the only ones for which the inverse function 
t = t(r) =B u ^(t) can be expressed in closed form. It is an algebraic function if 
a is any positive integer, and if a = 2 it can even be expressed in terms of radicals. 
But for any integer a > 1, its continuation to the r-plane is multiple- valued. 

No similar table of closed-form expressions seems to have appeared previously; 
though certain values of the traditionally normalized inverse incomplete beta 
function (especially for large (half-)integral values of a\ , c^ ) were calculated and 
tabulated long ago, owing to their importance in statistics ([Thompson! 1 1 94 lh . 
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